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Abstract. We have constructed equations of state involving various exotic forms 
of matter with large strangeness fraction such as hyperon matter, Bose-Einstein 
condensates of antikaons and strange quark matter. First order phase transitions 
from hadronic to antikaon condensed and quark matter are considered here. The 
hadronic phase is described by the relativistic field theoretical model. Later 
those equations of state are exploited to investigate models of uniformly rotating 
compact stars. The effect of rotation on the third family branch for the equation of 
state involving only antikaon condensates is investigated. We also discuss the back 
bending phenomenon due to a first order phase transition from K~ condensed to 
quark matter. 



1. Introduction 

Neutron stars are born in type II supernova explosions of massive stars 0. The matter 
in the neutron star interior is highly dense and cold. The matter density could exceed 
by a few times normal nuclear matter density Exotic components of matter such 
as hyperons, antikaon condensates and quarks may appear there. Such a dense and 
cold matter can not be found in laboratories. Therefore, neutron stars are very useful 
laboratories to study the properties of dense matter [3 E| ■ 

It is a challenging problem to constrain the equation of state (EoS) of dense matter 
in neutron stars. The theoretical investigation of mass-radius (M-R) relationship of 
compact stars is important because this can be directly compared with observations. 
Consequently the composition and EoS of dense matter may be probed. Though 
accurate masses have been measured for several neutron stars in radio pulsar binaries, 
we do not have enough informations about their radii. Observations of thermal x- 
ray and optical radiation from isolated neutron stars such as RX J185635-3754, may 
provide informations about their radii Another interesting candidate for radius 
measurement is x-ray bursters. Cottam et al. |S] estimated the gravitational redshift 
z = 0.35 of three spectral lines in x-ray bursts from EXO 0748-676. Recently, Villarreal 
and Strohmayer detected a 45 Hz rotational frequency for the same star 0. These 
observations lead to a radius of 9.5-15 km and mass l.h-2.'iM solar for EXO 0748-676 
|Hj . It is now being argued that the measurement of moment inertia of star A in double 
pulsar system PSR J0737-3039 where masses have been determined already, could give 
an accurate estimate of the radius and put important constrain on the properties of 
neutron star matter [7]. 
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Here we discuss various exotic forms of matter and their influence on equations 
of state and structures of static and rotating compact stars. In this context, 
we also report a new class of superdense stars called the third family in static 
IHl El El El Ej as well as rotating configurations JJj and the back bending 
phenomenon [2 13 El which could be a signature for a first order phase transition 
in the neutron star interior. In section 2, we describe the calculation of EoS and 
structure of compact stars. Results are discussed in section 3. And the outlook is 
presented in section 4. 

2. Model 

Here we describe models of static and rotating compact stars. Such models are 
constructed by solving Einstein's field equations, 

Gf"" ^ R'"' -^g''''R = 8TTT'"', (1) 

where R^'^ is the Ricci tensor, g^'^ is the space-time metric, R is the scalar curvature 
and T'^" is the energy-momentum tensor. The EoS is encoded in the energy- 
momentum tensor. The equilibrium structures of non-rotating stars are calculated 
from Tolman-Oppenheimer-Volkoff (TOV) equations. Here we are interested in the 
equilibrium structures of rotating compact stars. In this case, the metric has the form 

niin] 

ds^ = -e"'+Pdt^ + e^°'{dr^ + r^dO^) + e'^'Pr^sin^e{d(l) - Lodtf , (2) 

where metric functions depend on r and 0. Angular velocity of local inertial frames is 
denoted by uj. 

In the following paragraphs, we discuss two equations of state for dense matter. 
In the first case (case I), we consider a first order phase transition from hadronic 
matter including hyperons to K~ condensed matter followed by a second order 
condensation. In the other case (case II), we have first order phase transitions from 
nuclear to K~ condensed matter and then to quark matter. 

For case I, we have hadronic phase, antikaon condensed phase and a mixed of 
those two phases. The constituents of matter in the hadronic phase are all the species 
of the baryon octet, electrons and muons. Baryon-baryon interaction is mediated by 
the exchange of scalar and vector mesons. Here we adopt a relativistic field theoretical 
model to describe hadronic and antikaon condensed phase. The model is also extended 
to include hyperon-hyperon interaction through two additional strangeness meson 
/o(975) (hereafter denoted as a*) and (/)(1020). Therefore the Lagrangian density 
for baryon-baryon interaction is El E] 

Cb= '^^B (hfid^" - iriB + gaB<y - gujBlfj^UjP - ffpB7pts • P'') ■0-B 
B 

+ X! ida'BCr* - g^Bl,,<P^) V'B 
B 

- ^0^-0^" + \ml4>^4>'' ■ (3) 
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The scalar self-interaction term is 

The antikaon condensed phase is composed of baryons of the octet embedded in 
the condensate, electrons and muons. Similarly, we adopt a relativistic field theoretical 
approach for the description of (anti)kaon-baryon interaction. Here (anti)kaon- 
baryon interaction is treated in the same footing as baryon-baryon interaction. The 
Lagrangian density for (anti)kaons in the minimal coupling scheme is ^| ^] 

Ck = D*^KD''K ~ m*jlKK , (5) 

where the covariant derivative I?^ = 9^ + ig^K^^j, + ig<iiK4>ti + WpK^K ■ P^- The 
isospin doublet for kaons is denoted by K = {K^,K'^) and that for antikaons is 
K = (K~ , K'^). The effective mass of (anti)kaons in this minimal coupling scheme is 
given by 

m*ii ^ rriK - gaKO- - ga'Kcr* , (6) 

In the mean field approximation adopted here, the meson fields are replaced by their 
expectation values. The mean meson fields are denoted by a, cr*, wo, 4>o and po3- Each 
meson field equation in the antikaon condensed phase has contributions from antikaon 
condensates ^1^2- These equations differ from those of pure hadronic phase by the 
antikaon condensate terms. 

The in-medium energies of antikaon for s-wave (k ~ 0) is 

^K-,KO = m}^ - guK^o - 50X00 T ^ffpKPoa , (7) 

where the isospin projection /g^ = =Fl/2 for K~ meson and meson respectively. 
The onsets of antikaon condensations are given by [inilTT] . 

Mif- = Me , (8) 
^iKO = , (9) 

where nx- and //^o are respectively the chemical potentials of K~ and . 

The energy density (e'') and pressure (P^) in hadronic phase are related by 
= 'Yliit^i'^^i ~ ■ Similarly, the pressure in the antikaon condensed phase follows 
from the relation = ^ • pLiiii — , where rii is the number density of i-th species. 
The expressions for e'^ and are given as in Ref. 

The mixed phase of antikaon condensed matter and hadronic matter is governed 
by Gibbs phase rules and global conservation laws The Gibbs phase rules read, 
ph _ pK g^^^ ^ft _ .^^jjgj-g g^j.g chemical potentials of baryon 

B in hadronic and condensed matter respectively. The conditions for global 
charge neutrality and baryon number conservation are (1 — x)Q^ + xQ^ = and 
riB = (1 ~ x)"-B + X"-B where x is the volume fraction in the condensed phase. 
Charges in hadronic and K~ condensed phase are = qsn^ — n^^ — and 

= 'Y^g Qsn^ — uk- — Ue — respectively, where and rig are the number 
density of baryon B in pure hadronic and antikaon condensed phase and rij^- , and 
are number densities of antikaons, electrons and muons respectively. The total 
energy density in the mixed phase is given by e = (1 — + X^^ ■ 

For case II, the first order phase transition from nuclear matter to K~ condensed 
matter is treated in the same prescription as it is described above. However, the 
hadronic and antikaon condensed phases are composed of neutrons, protons, electrons 
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Figure 1. The equation of state, pressure P versus energy density e, for two 
cases as described in the text. 



and muons. In this case, we also consider a first order phase transition from antikaon 
condensed matter to quark matter. The pure quark matter is composed of u, d and s 
quarks. We describe the EoS of pure quark matter in the MIT bag model |20|. The 
free energy density of non-interacting quarks is given by, 

" = ^ E f''^dpp\^p^ + m'i-fi,) + B, (10) 

i—u.d,s ^ 

where B is the bag energy density. The pressure in the quark phase is given by 
pQ = -fl and the energy density (e*^) is obtained from the Gibbs-Duhem relation. 
And the mixed phase is governed by the Gibbs phase rules and conservation laws. 

3. Results &: Discussions 

In this calculation, we adopt GMl parameter set pi where nucleon-meson coupling 
constants are determined from the nuclear matter saturation properties. The vector 
meson coupling constants for (anti)kaons and hyperons are determined from the quark 
model ^Hl whereas the scalar meson coupling constants for hyperons and antikaons 
are obtained from the potential depths of hyperons and antikaons in normal nuclear 
matter ^2^]. As the phenomenological fit to the K~ atomic data yielded a very 
strong real part of antikaon potential U/^ = —180 ± 20 MeV we perform this 
calculation with an antikaon optical potential of -160 MeV at normal nuclear matter 
density (no = 0.153/m^'^). The coupling constants for strange mesons with hyperons 
and (anti)kaons are taken from Ref.^JEI- Also, we consider a bag constant B^^'^ = 
200 MeV and strange quark mass = 150 MeV for quark EoS in case II. 
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Figure 2. Mass-radius relationship of non-rotating stars calculated with the 
equation of state for case I. 



The equations of state, pressure versus energy density, for case I and case II are 
exhibited in Figure 1 by soUd and dashed hues respectively. For both cases, the upper 
and lower boundaries of the mixed phase of hadronic and K~ condensed matter are 
identified by two kinks. The mixed phase begins at energy density e — 6.11x10^^ 
g/cm^ (or 2.23no). For case I, A hyperons are populated in the mixed phase at 2.51no 
whereas S~ and hyperons are populated at much higher densities. A second 
order condensation occurs in case I just after the mixed phase is over at e = 
1.17x10^^ g/cTP?. For case II, the phase transition to quark matter is delayed due to 
the early onset of K~ condensation depending on the parameters. The first order K~ 
condensed-quark matter phase transition begins dX e — 1.59x10^'^ g/cw? and ends at 
e = 2.77x10^^ g/cTo?. In this case, the extent of the mixed phase is quite large. We 
note that the overall EoS for case I is softer than that of case II. 

In Figure 2, the mass-radius relationship of non- rotating stars is displayed for 
the EoS of case I. The filled circles correspond to the maximum masses. We obtain 
a maximum neutron star mass of l.blM solar corresponding to a radius of 12.8 km 
as shown in the figure. It is found that after the positive slope neutron star branch, 
there is an unstable region followed by another positive slope compact star branch. 
This new branch of compact stars is the result of the discontinuity in the speed of 
sound due to the kinks or first order phase transition in the EoS for case I. This new 
sequence of superdense stars beyond the neutron star branch is called the third family 
branch. From the study of fundamental mode of radial vibration, it was found that 
the third family branch was a stable one |11| . The compact stars in the third family 
branch have different compositions and smaller radii than their counterparts on the 
neutron star branch 11, .13: . The maximum mass of the compact star in the third 
family branch is l.bbMsoiar corresponding to a radius 10.7 km. 
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Figure 3. Gravitational mass as a function of equatorial radius for case I. 
The static and mass shedding limit sequences along with fixed angular velocity 
sequences are displayed here. Also, normal and supramassive sequences are 
plotted. 



Now we discuss the mass-equatorial radius relationship of rotating sequences 
calculated with the EoS of case I. This is shown in Figure 3. In the preceding 
paragraph, we have discussed the third family of compact stars in the static sequence 
calculated with this EoS. It is interesting to watch what happens to this third family 
branch in rotating cases. In Figure 3, the extreme left and right curves denote the 
static and mass shedding limit sequences. Each compact star on the mass shedding 
limit rotates with a maximum frequency called the Kepler frequency and it becomes 
unstable beyond this frequency. Fixed angular velocity curves are also shown in this 
figure. On each curve, the maximum mass in the neutron star branch is denoted by 
solid circles whereas the maximum mass in the third family branch is represented by 
open circles. We observe that with increasing rotation, the maximum star masses also 
increase and are shifted to smaller central energy densities or larger radii. This is 
attributed to the presence of centrifugal force in rotating stars. On the other hand, 
we note that the third family of compact stars disappears beyond the fixed angular 
velocity curve O = 2158 sec~^. The reason is the strengthening of the centrifugal force 
with increasing rotation. Consequently, fast rotating stars could not probe the high 
density part of the EoS which is responsible for the origin of a third family branch. 

In Figure 3, we also plot fixed baryon number evolutionary sequences. Here two 
normal sequences are shown by curves I and II and two supramassive sequences by 
curves III and IV. Fixed baryon number sequences provide important informations 
about isolated rotating stars as they slow down with time due to the loss of energy 
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Figure 4. Moment of inertia (I) is plotted with angular velocity (n) for the 
EoS of case II. Normal and supramassive sequences are shown by curves I and II 
respectively. The other curve is the mass shedding limit sequence. 



and angular momentum through electromagnetic and gravitational radiation. We 
observe from the figure that non-rotating stars are the members of normal sequences. 
However, no member of supramassive sequences is a non-rotating compact star. When 
a supramassive star loses sufficient angular momentum such that the centrifugal 
support against self gravity is impossible, it would then collapse to a black hole. 
This was considered as a possible scenario for gamma ray bursts pS] . 

We also investigate the mass-equatorial radius relationship of the static, mass 
shedding limit and fixed angular velocity sequences calculated with the EoS for case 
II. However, we do not find any third family solution in the static and rotating cases. 
But we find some interesting results from the study of a supramassive sequence. We 
display moment of inertia (I) with angular velocity (f2) in Figure 4. Curve I implies 
the normal sequence with baryon number iV^ = (2.15 x 10^^) whereas curve II is a 
supramassive sequence with Nf, = (2.45 x 10^^). The other curve is the mass shedding 
limit sequence. For the normal sequence, the moment of inertia always decreases 
with decreasing angular velocity. However, curve II exhibits that after the initial spin 
down of the neutron star along this sequence, there is a spin up followed by another 
spin down. This is known as the back bending phenomenon. We observe that this 
back bending phenomenon is the result of the first order phase transition from antikaon 
condensed to quark matter. It is worth mentioning here that the supramassive neutron 
star is stable on the back bending segment jj^ ■ We also note that the period of the 
neutron star in the back bending part varies from 1.10 ms to 1.35 ms corresponding to 
the variation in angular velocity from 5700 s~^ to 4660 s~^. The minimum observed 
pulsar period is so far 1.56 ms. Future pulsar surveys would tell us whether pulsars 
rotating faster than the currently known one exist in nature or not. 
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4. Outlook 

We have investigated equations of state including exotic forms of matter and computed 
the structures of static and rotating compact stars. Many interesting results such as 
the third family of compact stars containing exotic matter, the back bending effect 
which could signal a first order phase transition in the neutron star interior, are 
obtained here. Our results are sensitive to input parameters such as antikaon optical 
potential depth and bag constant. Future experiments on deeply bound states of 
antibaryon-nucleus at GSI, Germany and _R^-nucleus at Japan Proton Accelerator 
Complex (J-PARC) will provide opportunity to investigate dense and cold matter 
in laboratories and determine those parameters accurately. It would be possible to 
constrain the composition and EoS of dense and cold matter in compact stars from 
observations and experiments in future. 
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